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We investigate the kinetics of nonlinear collision-induced fragmentation. We obtain the fragment 
mass distribution analytically by utilizing its travelling wave behavior. The system undergoes a 
shattering transition in which a finite fraction of the mass is lost to infinitesimal fragments (dust). 
The nature of the shattering transition depends on the fragmentation process. When the larger of 
the two colliding fragments splits, the transition is discontinuous and the entire mass is transformed 
into dust at the transition point. When the smaller fragment splits, the transition is continuous with 
the dust gaining mass steadily on the account of the fragments. At the transition point, the fragment 
mass distribution diverges algebraically for small masses, c(m) ~ m~", with a = 1.20191 . . . 
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I. INTRODUCTION 

Fragmentation occurs in numerous physical phenom- 
ena and industrial processes Ems, 11,13 Examples 
include breakup of liquid droplets [g and atomic nuclei 
[t'I, polymer degradation shattering of solid objects 
|9, 10], meteor impacts, and mineral grinding. Idealized 
models of such physical phenomena are also useful con- 
ceptual tools for describing complex systems such as fluid 
turbulence, spin glasses geneticpopulations [l^lTsj . 
and random Boolean networks 

In some cases, for example in polymer degradation, the 
evolution of a fragment depends only on its size. There- 
fore, fragments do not interact and such processes are 
inherently linear. In other cases including grinding pro- 
cesses, explosions in an enclosed volume, and breakup of 
eddies in a turbulent flow , interactions between frag- 
ments are essential. Such fragmentation processes are 
intrinsically nonlinear |17l 18, 19, 20] . In this study, we 
show that the nature of the mass distribution changes 
qualitatively due to nonlinearities. 

We investigate a basic class of nonlinear fragmentation 
processes where binary collisions are the cause of break- 
age. We show that such processes exhibit a shattering 
transition where infinitesimal fragments (dust) carry a 
finite fraction of the mass in the system. We consider the 
simplest realization where one of the two colliding frag- 
ments breaks into two pieces. Generically, the number of 
fragments diverges in a finite time, indicating shattering 
into dust. 

The nature of the shattering transition depends sensi- 
tively upon the details of the fragmentation process, in 
particular, which of the two colliding particles splits. We 
investigate three possibilities: (A) either, (B) the larger. 
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and (C) the smaller of the two particles breaks into two 
fragments upon collision. In the first two models, as the 
transition occurs the entire mass is instantly transformed 
into dust. In the third model, the dust mass gradually 
increases once the shattering transition occurred. 

In contrast with linear fragmentation processes, ex- 
plicit solutions of the nonlinear and non-local rate equa- 
tions are generally not possible. Nevertheless, the most 
important physical characteristics can still be obtained 
analytically. Interestingly, the fragment mass distribu- 
tion attains a travelling wave form as the transition is 
approached. Of the spectrum of possible propagation ve- 
locities, the extremal one is selected and it characterizes 
typical and extremal behaviors of the mass distribution. 
In the case of model C, at the shattering transition, the 
mass distribution is algebraic for small masses, with a 
transcendental exponent. Past the transition, the frag- 
ment mass distribution approaches a universal form. 

We first consider the number density that manifests 
the shattering transition (section II). Then, we analyze 
the fragment mass distribution using rate equations for 
a deterministic version (section III) and a stochastic ver- 
sion (section IV) of the fragmentation process. Finally, 
we summarize our results and outline a few suggestions 
for future work (section V). 



II. THE NUMBER DENSITY 

Consider a fragmentation process where at each (bi- 
nary) collision event, one particle splits into two pieces 
while the second particle remains intact. Without loss of 
generality, the collision rate is set to unity. Analogous to 
the kinetic theory description of collisions in molecular 
gases, we assume perfect mixing, namely, absence of spa- 
tial correlations between fragments. The total fragment 
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density, N{t), evolves according to the rate equation 

±N{t)^N\t). (1) 

Without loss of generality, the initial density is set to 
unity, A^(0) = 1, and therefore, the total density is 



Nit) = 



1 



1 -t 



(2) 



In a finite time, the number of fragments diverges and 
the average fragment mass vanishes. This divergence in- 
dicates that the system undergoes a shattering transition 
a.t tc — I. 

Let 7" = /q dt' N{t') be the average number of collisions 
experienced by a fragment up to time t. This quantity 
diverges logarithmically 



IniV(t) = In 



1 



1 - t 



(3) 



This "collision counter" provides a convenient alternative 
measure of time. 



III. DETERMINISTIC FRAGMENTATION 

To complete the model definition we have to specify 
which of the fragments splits, and how it splits. Follow- 
ing Cheng and Redner [l^, we consider three possibili- 
ties: (A) a randomly chosen, (B) the larger, and (C) the 
smaller fragment splits upon collision. In this section, 
we consider a deterministic rule where fragments split 
into two equal pieces. In the next section, we show that 
stochastic rules result in qualitatively similar behaviors. 



A. Random particle splits 

We start with the case where a randomly selected par- 
ticle splits upon collision (this is equivalent to having 
both particles split). For simplicity, we focus on monodis- 
perse initial conditions where all particles have unit mass, 
TO = 1. Then, a fragment produced by n collision events 
has mass to = 2~". Let c„(t) be the density of such 
fragments at time t. This density evolves according to 



-c„(t) = 7V(t)[2c„_i(<) 
at 



with the total density N{t) — '^'jLoCjit)- Summing up 
Eqs. Q we indeed recover Eq. (^). Also, the total mass, 
M{t) = X;^o Cj{t), is conserved, M{t) = 1. 

In terms of the collision counter, the process is linear, 
^ c„ = 2c„_i — c„, and subject to the monodisperse 
initial conditions Cn(0) — Snfij the exact solution is the 
Poissonian density |18| 



c„(t) 



(2t)" 



(5) 



At the shattering time tc = 1 (corresponding to r = oo), 
the densities vanish: c„(i = 1) = for all n. There- 
fore, the fragment mass density undergoes a first-order 
(discontinuous) transition, M{t) = 8(ic — t) with 8 the 
Heaviside step function. In other words, the entire mass 
is shattered into dust and there are no particles with pos- 
itive mass mil EH. 

Near the shattering transition, i.e., as r — > oo, the 
mass distribution approaches 



c„(t) 



VT 



(6) 



where v = 2 and G{x) = (27r)"i/2 exp{-x^/2) is the 
Gaussian distribution. Since n — log2(l/TO), the mass 
distribution becomes log-normal, a behavior typical to 
fragmentation and cascade processes 0, 0, [T^ |21| . 



B. Larger particle splits 

Now in a collision the larger particle splits into two 
equal pieces. If the colliding particles have the same 
mass, a randomly chosen particle splits. The fragment 
mass density, c„(i), satisfies the rate equation 



at 



2c^, 



- 



(7) 



where An is the cumulative density of fragments of mass 



j{t). The initial condi- 



2 " and smaller, An(t) 
tions are c„(0) = (Sn,o- One can verify that the mass is 
conserved, M{t) — 1, and that the total density is given 
by Eq. ©. 

The density co(t) of unit mass particles satisfies the 
Bernoulli equation, ^ cq = Cq — 2coN. Using Eq. Q and 
the initial condition co(0) = 1 gives 



co{t) = 



3(1 -t)^ 
2 + {l-tf 



(8) 



For sufficiently small n, one can obtain the leading 
asymptotic behavior near the shattering transition. Since 
An N a.s t ^ 1 and the last two terms on the right- 
hand side of Eq. are asymptotically negligible, the 

- c„) which 
Therefore, 

(9) 



rate equations simphfy to ^ c„ = 2iV(2c„_i - 
(4) are identical (up to the factor 2) to Eqs. 



c„(t) oc e' 



-2. W 



Apart from logarithmic corrections, the densities vanish 
quadratically: c„(i) oc (1 — i)^. We conclude that the 
shattering transition remains discontinuous (see Fig. 
Figure suggests studying the normalized distribution 
N~^Cn{t). Below, we show that as a function of t, the 
normalized fragment mass distribution follows a universal 
behavior in the large-n limit. 



3 





FIG. 1: The normalized fragment size distribution. Shown is 
N^^Cn{t) versus t for n — 0, 1, 2, 4, and 6. 



FIG. 2: The traveUing wave. Shown are numerical solutions 
of Eq. (|nj for n = 2, 4, 6, and 8. 



The rate equations simplify in terms of the cumu- 
lative densities: 



2Ai 



(10) 



This equation holds for = ^ if we set A^i = Aq. 
The initial conditions are ^n(O) = (^n,o- We characterize 
time by the collision counter Q and normalize the size 
density by the total number density, F„(r) ~ N^^An{t). 
These transformations yield 



(11) 



Asymptotically, this equation admits a travelling wave 
solution Fnir) /{n — vr) as shown in Fig.|21 The wave 
form f{x) satisfies the difference-differential equation 



ax 



f{x)+f{x)-2f\x-l), (12) 



and is subject to the boundary conditions f{—oo) = 1 
and /(oo) = 0. Remarkably, the velocity v can be deter- 
mined without solving the nonlinear and non-local dif- 
ferential equation (|12|) exactly. It follows from the expo- 
nential behavior attained by f{x) far behind the front: 
1 — f{x) ^ e^^ as X — oo. Together with Eq. ifT^ it 
yields a "dispersion" relation between the velocity v and 
the decay coefficient A, 



3-4:6' 



(13) 



Out of the spectra of possible velocities v e (— oo, Wmax]j 
the maximal value is selected. At the maximum, we 
have 3e-^ = 4(1 + A), from which A = 0.961279 and 



V = 1.52961. Alternatively, the velocity is the smaller 
root of wln(4e/?;) = 3. 

Velocity selection underlies numerous situations, yet it 
has been rigorously established only for a few non-linear 
parabolic partial differential equations, t ypically occur- 
ring in reaction-diffusion problems |2a . l2al24Ll2al2al27l| . 
Recently, velocity selection has been also applied to a 
host of difference and difference-differential equations 
'28', "2^, "3^, I31L [3^ including a linear fragmentation pro- 
cess ^0} . Typically, the selected velocity gives key phys- 
ical characteristics such as t he g rowth velocity of a sur- 
face in deposition processes the extremal heights 
of random trees 32J. 

The typical behavior of the fragment mass density fol- 
lows from the travelling wave form 



Ng{n 



(14) 



with g{x) = f{x) — f{x + 1). The front location w 
VT characterizes typical fragments and the typical mass 
TO* — 2~"* shrinks as 



(i-tr 



(15) 



with a — vln2 ^ 1.06024 as i ^ 1. The typical mass de- 
cays slower than in model A where ct = 21n2 = 1.38629. 
Another difference between models A and B is manifested 
by the width: In contrast with the diffusive broadening 
in model A the width saturates at a finite value in model 
B. Yet, fundamentally the shattering transitions are the 
same in both models — the entire system is instantly 
transformed into dust at the transition point. 

The extremal behavior of the fragment mass density 
follows from the tails of f{x). The behavior far ahead 
of the wave front {x 00) is a sharp double exponen- 
tial decay, as implied by the leading terms in Eq. H12|l . 
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v^f{x) = —2/^ (a; — 1). In summary, the extremal be- 
haviors are 



1 - Ci e-^^ 

2^ exp(-C2 2^) 



X 
X 



-oo, 

-CO. 



(16) 



We now re-express the mass distribution in terms of the 
ordinary mass variable m = 2^". The two distributions 
are related via c{m)dm = Cndn. Near the shattering 
transition, the mass distribution attains the scaling form 

c(to) -^jr (^-HL^ . Equation (O leads to the following 

extremal behaviors of the scaling function: 



^(z) 



■ exp I 



2 Z 



Z > 1, 

z < 1; 



(17) 



with a — 1 -|- A/ In 2 = 2.38683. Hence, large masses 
(relative to the typical mass) are suppressed algebraically, 
while small masses are suppressed exponentially. 

Generally, in fragmentation processes the mass distri- 
bution has a scaling form and this is indeed the case for 
collision-induced fragmentation. However, the nonlinear 
nature of the process results in qualitative changes to the 
scaling behavior. The similarity solutions have two scales 
characterizing the front location and fluctuations around 
it in the linear case (model A). In contrast, only a single 
scale underlies similarity solutions in the nonlinear case 
(model B). 



C. Smaller particle splits 

When the smaller particle splits upon collision the frag- 
ment size densities satisfy the rate equations 



dt 



(18) 



where Bn = X]J=o ^he cumulative density of parti- 
cles with mass larger than 2~". 

The density of unit mass particles is readily found by 
solving Co — "Cg. The next density can be found as well 



co(i) = 



1 



1 + t ' 
2 {l + tf 



(19) 



l + t 2{l + tY + l 



These explicit results already demonstrate that densities 
are positive at all times. Hence, the total mass density 
Mit) also remains positive after the shattering transition. 

The kinetics just below and at the shattering transition 
can be determined using the travelling wave behavior. 
The cumulative distribution obeys ^ i?„ — 2Bf^_^ — 
which is identical to Eq. (|1U|) : the initial conditions, how- 
ever, are different: Bn{0) = 1 — (5n,o- The transformed 
distribution F„(t) = N'^B^ again evolves according to 
Eq. (|ll|l . Asymptotically, it admits a travelling wave so- 
lution, Fnir) — > f{n — vt), with the wave form f{x) sat- 
isfying Eq. H12|) . However, the boundary conditions are 




m 



FIG. 3: The mass distribution at the shattering time. Nu- 
merical integration of the rate equations 11811 are compared 
with the theoretical prediction H23|l . 



reversed, F{—oo) = and F{oo) = 1, leading to different 
quantitative and qualitative results. 

Both extremal behaviors are now exponential 



1 



X 
X 



-OO, 

-oo. 



(20) 



The behavior far ahead of the front is used to determine 
the velocity. The dispersion relation is 



4p/ - 3 
A 



(21) 



and the extremum selection principle gives A = 0.58013 
and V ^ 7.14509. Numerically, we confirmed this velocity 
to within 0.01%. Interestingly, v is the larger root of the 
same (as in model B) equation vln{Ae/v) = 3. We note 
that the velocities satisfy vb < va < vc- 

The fragment size distribution follows the travelling 
wave form H14|l with 17(2:) ~ f{x + 1) — f{x). The typical 
mass shrinks according to H15|l with cr = i>ln2 = 4.9526 
near the shattering point. The exponential tails of the 
wave form imply algebraic tails for the scaling function 
underlying the mass distribution 



z>l, 

-r''^ z < 1; 



(22) 



with a = 1 + (wln2)-i 9i 1.20191 and /? = 1 - A/ In 2 = 
0.163049. 

Our major result is that the mass distribution diverges 
algebraically at the transition time [33l |: 



c(m, 1) ^ TO 



(23) 



for TO — > with the transcendental exponent a = 1.20191 
(Fig. 3). This behavior can be obtained from the large-z 
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FIG. 4: Fragment versus dust mass. Shown are the fragment 
mass M{t) (solid hne) and the dust mass fi{t) (dashed hne) 
versus time t. 



FIG. 5: The asymptotic behavior of the size density. Shown 
is 2-"/^Cn(i) versus t ior n = 6, 10, 20, and 30. 



behavior of ^(z) . Although in general the travelling wave 
form implies time-dependent densities, when z — > oo, the 
mass densities become stationary. 

Model C exhibits a rich post-transition behavior. The 
explicit solutions ifT^ suggest that c„ ~ 7„ when t 
oo. Indeed, this behavior is compatible with Eqs. H18|) 
and the cumulative amplitudes r„ = J2^=o satisfy the 
recursion relation - r„ = 2T'^^_^ with Tq = 1. The 
amplitudes grow exponentially, 7„ ^ r„ ^ 2"/^. Sum- 
ming over densities, the total fragment mass decays as 

M{t)~Ct-^ as t^oo, (24) 

with C = Er=o2""7n = 2.66084. Thus, the total 
fragment mass remains positive at all times. The dust 
mass, ii{t) = 1 — M{t), vanishes at the shattering time, 
/i(l) = 0, and it gradually increases for f > 1 (Fig. 0J. 
Only in the long time limit it accounts for the entire 
mass in the system. We conclude that in model C, the 
shattering transition is continuous. 

Numerically, we observe that for sufficiently large n, 
the densities follow a universal behavior (Fig. 0) 

Cn{t)^2"/^u{t). (25) 

While this ansatz is asymptotic with respect to n, it holds 
for all times. The function u{t) vanishes below the shat- 
tering time and grows linearly afterwords, u(t) ~ (i — 1) 
for t — 1 ^ 0. Hence, this function plays the role of an 
order parameter. Note also that u{t) ~ t~^ as i ^ oo. 

The order parameter and the total dust mass are inti- 
mately related. Consider the total mass density of frag- 
ments of mass 2"*^ or larger: M^^^(t) = X]n=o ''"(^)- 
From the rate equations H18|) , this mass density decreases 



according to 

^M«(i) = -2-'=Cfc(2i?fc + Cfc). (26) 
at 

The flux of mass from fragments into dust is simply 
^ = - hmfe_oo jrt Using Eq. ^I^, the right-hand 

side of Eq. (|26|l approaches (3 + 2\/2 ) M^(i) in the limit 
k oo and therefore, 

^ti{t) = {3 + 2V2)u^it). (27) 

This in turn shows that the dust mass grows according 
to fi{t) {t — 1)^ past the transition. 



IV. STOCHASTIC FRAGMENTATION 

We now briefly describe a generalized collision-induced 
fragmentation process where splitting is stochastic. 
Specifically, a particle of mass m splits into two frag- 
ments of mass m' and m — m' with m' chosen stochas- 
tically from the interval < m' < m according to some 
fixed distribution. We focus on the simplest case of uni- 
form splitting, i.e., m' is chosen uniformly in [0,to]. 



A. Model A 

When a randomly selected particle splits, the mass 
density c(m, t) satisfies 

— c(m,T) = -c(to,t) + 2 / — c(m',T). (28) 
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The kernel 1/m' reflects the uniform splitting probabil- 
ity and the collision rate N is absorbed by the colli- 
sion counter t. This equation is solved using the Mellin 
transform and for the monodisperse initial condition, 
c(m, 0) — 5{m — 1), one finds 



c{m,T) = e ^(5(m— l) + e 



/ 2t 
h^ 



h 



Srln — 
m 



with Ii the modified Bessel function. The first term on 
the right-hand side simply describes the density of par- 
ticles that have yet to collide. The second term simpli- 
fies asymptotically. Making the transformation m = e^" 
leads to a normal distribution as in © with the propa- 
gation velocity v = Thus, the qualitative behavior is 
independent of the splitting rule. 



B. Model B 

When the larger of the two fragments splits, the rate 
equations for the mass density are 

d f°° dm' 

■^c{m) = 4 / — -i-c{m')A{m') - 2c{m)A{m) , (29) 

with the cumulative density A(m) = dm'c{m'). We 
employ the same transformations used in the determin- 
istic case. Characterizing the mass m by "index" n via 
m — e~", the fragment size density, c(n), evolves accord- 
ing to 

c{n) = 4 / dn'e" ^''c{n')A{n') ~ 2c{n)A{n) (30) 



dt 



with A{n) — dn'c{n'). This cumulative distribution 
satisfies 

^^A{n) = A\n) - 2 dn' ~n^A\n') . (31) 

Expressing time in units of the collision counter and nor- 
malizing by the total density, F{n,T) = N~^A{n), we 
transform Eq. (|31|) into 



Fin) ^ F^{n) ~ F(n) - 2 / dn' e" -" ^^^(n') , 

Jo or' 



Seeking a travelling wave solution F(n,T) f(n — vt) 
yields the non-linear integro-differential equation 

v^f{x)^f{x)-f{x) + 2 £^dyey-^^f{y), (32) 

subject to the boundary conditions /(— oo) ~ and 
f{oo) = 1. The exponential decay 1 — f{x) ^ exp(Ax) as 
X — > oo gives the dispersion relation v — 4(1-|-A)^^ — A^^ 
and the extremum selection principle yields A = u = 1. 
Close to the shattering transition, the typical mass is pro- 
portional to the average mass, m, ~ (1 — t) [s^. The 
mass densities behave as in the deterministic case and 
the extremal behaviors H17|l are recovered with a — 2. 
The nature of the transition is discontinuous, as in the 
deterministic case. 



C. Model C 

When the smaller particle splits upon collision, the rate 
equations for the mass density are 

3 dm' 

^c(to) = 4 / c{m')B{m') - 2c{m')B{m') (33) 

dt to' 

with B{rn,t) — dm' c{m' , t) . In terms of the index 
n, the cumulative density B{n,t) = dn'c{n') obeys 
Eq. (|31|l . The normalized cumulative density again ad- 
mits the travelling wave form. The velocity and de- 
cay rate are v — 9 and A = i. At the shattering 
time, the (finite) mass distribution diverges algebraically: 
c(to, 1) ^ TO~" with a = 10/9. Past the shattering 
transition, the asymptotic ansatz c(m, t) —^ m-^/'^uit) 
holds for small masses and the dust mass is related to 
the order parameter via fJ- = . We conclude that 
qualitatively, the shattering transition is similar to the 
deterministic case. 



V. DISCUSSION 

We investigated kinetic properties of collision-induced 
fragmentation processes. Generally, the mass is trans- 
ferred from finite fragments into infinitesimal dust in a 
finite time. The nature of the shattering transition de- 
pends on the fragmentation process. When the larger of 
the colliding particles splits or when a randomly selected 
one splits, the transition is discontinuous and the entire 
mass is transformed into dust instantaneously. When 
the smaller particle splits, the transition is continuous, 
with the dust accumulating gradually past the shatter- 
ing transition. In this case, finite fragments always carry 
a non-zero fraction of the mass. 

Model A is essentially linear and thus solvable. For 
models B and C the nonlinear and non-local governing 
equations can not be solved in a closed form. Neverthe- 
less, in the vicinity of the shattering transition we were 
able to obtain the most important characteristics analyt- 
ically by utilizing the traveling wave form of the fragment 
mass density. The mass distribution follows a scaling be- 
havior with a single characteristic scale, in contrast with 
the two scales found for linear processes. 

For model C, the post-shattering behavior is nontriv- 
ial. At the transition point, the mass distribution decays 
algebraically, c(to) ^ m^", with a transcendental expo- 
nent a = 1.20191 ... in deterministic fragmentation and 
a rational exponent a = 10/9 in stochastic fragmenta- 
tion. We have also demonstrated that the mass densities 
exhibit universal asymptotic behavior (|25|l in the post- 
shattering region. 

A challenging open problem is the complete post- 
shattering behavior in model C, for example, the time de- 
pendent dust mass. This is largely a mathematical prob- 
lem since physically, the breakage of sufficiently small 
fragments is impossible. For instance, micro-cracks on 
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the surface of the fragment are often precursors for break- 
age. The number of such surface defects is proportional 
to the surface area, so sufficiently small fragments are 
effectively unbreakable. 

We focused on the leading asymptotic behavior. There 
are however corrections to the linear front propagation 
[H, m, HI, [231 . The traveling wave solution is actu- 
ally a function of x = n — X{t) with the position of the 
front X(r) given by X{t) =vT±^\nT + 0{1). The 
plus and minus signs correspond to model B and C, re- 
spectively. This translates to a logarithmic correction to 
the typical mass ((T3|) . 

We treated the problem using a mean-field rate equa- 
tion approach. Thus we ignored correlations between the 
colliding particles. In principle, spatial correlations may 
be important up to some critical dimension beyond which 
they can indeed be ignored. The analysis of this possibil- 
ity requires a more complete description of the process. 
Particularly, one must specify the transport mechanism. 

Collision-induced fragmentation arises most naturally 
in processes where particles moves ballistically between 



collisions. Using dimensional analysis we argue that the 
shattering transition always occur in ballistic fragmenta- 
tion. The typical mean free time T, velocity particle 
cross section s, and number density N are related via 
NvTs ~ 1. Mass conservation implies m ~ N~-^ (here 
TO is the typical mass), while energy conservation gives 
1; ~ 1. Finally m - s''/(''-i) yields s - N'^+^Z'^. In 
particular, T ^ N"^/'^. The particle density evolves ac- 
cording to ^ iV = N/T, or ^ iV - N^+i/d from which 
N ^ {tc — t)~'^ . This suggests that in ballistic fragmen- 
tation, the shattering transition occurs in arbitrary di- 
mension d. Using effective d-dimensional collision rates 
((X N^^'^) one can convert the "one-dimensional" results 
in this study into a general mean-field theory. 
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